Model-Independent Reionization Observables in the CMB 
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We represent the reionization history of the universe as a free function in redshift and study the 
potential for its extraction from CMB polarization spectra. From a principal component analysis, 
we show that the ionization history information is contained in 5 modes, resembling low-order 
Fourier modes in redshift space. The amplitude of these modes represent a compact description 
of the observable properties of reionization in the CMB, easily predicted given a model for the 
ionization fraction. Measurement of these modes can ultimately constrain the total optical depth, 
or equivalently the initial amplitude of fluctuations to the 1% level regardless of the true model for 
reionization. 



I. INTRODUCTION 

The WMAP experiment has recently detected the 
reionization of the universe through the large-angle po- 
larization of the cosmic microwave background (CMB). 
Rescattering of CMB radiation bearing a quadrupole 
temperature anisotropy leads to a small linear polar- 
ization which is then correlated with the temperature 
anisotropy itself. It was through this cross-correlation 
that WMAP measured the total optical depth to Thom- 
son scattering, r = 0.17 ± 0.04 0. 

The CMB power spectra potentially contain more in- 
formation than the optical depth integrated over the 
whole ionization history 0. Unfortunately, the nature 
of the ionizing sources is not currently well-understood 
(e.g. H|, c.f. |4|), and therefore the detailed ionization 
history should not be treated in the traditional CMB ap- 
proach of adding well-motivated model parameters to a 
parameter estimation chain. Attempts to do so can re- 
sult in biases and conflicting results. For example, using 
an ionization history derived from a particular numerical 
prescription [{| it was found that the partially ionized 
epoch was undetectable in the CMB, while analyses us- 
ing semi-analytic models 0, showed that the predicted 
multiple epochs of reionization and partial recombination 
can be detected 0,0. As another example, by assuming 
an overly simplistic step function model for reionization, 
the total optical depth r in the semi-analytic models can 
be mis-estimated by up to a few 10 -2 [!j; while not a 
large amount it is already approaching the statistical er- 
rors in the first year WMAP data. 

These uncertainties argue for a more model- 
independent approach to the phenomenology of reioniza- 
tion in the CMB. In this paper, we begin by considering 
a complete basis for the ionization history in fjn] In AIIII 
we employ a principal component analysis to isolate the 
CMB reionization observables. Throughout we work in 
a fiat ACDM cosmology with h = 0.72, Q m h 2 = 0.145, 
f2f,/i 2 = 0.024 and n = 1. In the absence of system- 
atic errors, most of the ionization history information 
comes from the power spectrum of the E-mode polar- 
ization Cf E as opposed to the temperature-polarization 



cross correlation measured by WMAP [l(| and for sim- 
plicity we will here consider only this power spectrum. 
Nevertheless the main results are applicable to the cross 
power spectrum as well. 



II. MODEL-INDEPENDENCE 

To zeroth order, the total Thomson optical depth pa- 
rameterizes the effects of reionization independently of 
the specific model for the ionization history. Firstly, the 
acoustic peaks in the CMB power spectra are lowered 
by e~ T in amplitude. Because the effect is degenerate 
with a change in the initial normalization of the curva- 
ture fluctuations 5^ , it is absorbed into fixing the combi- 
nation 8qe~ t . Secondly, rescattering generates large an- 
gle temperature fluctuations through the Doppler effect 
and thirdly, it generates large angle polarization. Both 
effects have an rms that depends mainly on r but since 
the primordial polarization is free of large angle contribu- 
tions, it is the better probe of reionization. Finally there 
are further signals from inhomogeneous scattering that 
appear beyond the damping tail (£ > 2000) but whose 
phenomenology is inextricably tied to the model for the 
ionizing sources; we will not consider them further here. 

In detail, the residual temperature and polarization ef- 
fects depend on the reionization history, mainly through 
a sensitivity to the horizon scale at the epoch of rescat- 
tering |ll| . We assume the hydrogen ionization fraction x 
to be a free function of redshift. A complete basis for this 
function can be formed from a series of delta functions in 
redshift. We approximate this series as Kronecker delta 
functions 5x(zi) spaced by Az with linear interpolation 
between the points (see Fig.QJ. The spacing is chosen to 
capture all of the information in the CMB and also be 
sufficiently fine to reproduce features in the ionization 
model. We take Az = 1 throughout. 

A specific model can then be represented as a sum 
over the delta function modes in the approximation of 
linear interpolation between the points. Consider a fidu- 
cial ionization history with a hydrogen ionization fraction 
of x = 1, Zi < 7; x = 0.5, 7 < Zi < 26; x — x ICC , Zi > 26 
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where a; roc (w 0) is the ionization fraction coming out of 
standard recombination (see Fig. 0. Compare that to 
a step function ionization with x = 1, z\ < 17 and x lec 
Zi > 17 (Fig.[IJ. Both models have r « 0.17 but differ in 
the coherence scale of the polarization. These two models 
are clearly distinguishable in principle (see Fig. top) 
and so more information than the optical depth can be 
extracted from precise measurements 0, . 

The parameters x(zi) may be simply appended to the 
usual CMB parameter estimation chain but that ap- 
proach has several drawbacks. The Az spacing of the 
points is arbitrary and moreover for Az = 1 the param- 
eters are highly degenerate and would cause numerical 
problems in a likelihood analysis. To see this, consider 
the effect of a delta-function perturbation to the fiducial 
ionization history given above and quantify the response 
in the power spectrum as a transfer matrix 



d\nCf E 

dx(zi) 



(1) 



where we vary 7 < 2, < 26 as motivated by the range 
in redshift over which reionization is expected to occur 
[3- The derivative is calculated by a double-sided fi- 
nite difference of Sx = 0.05. This matrix may be viewed 
as a transfer function for small perturbations from the 
fiducial model. The main feature is that perturbations 
at relatively low redshift appear at low £ and those at 
high redshift at high I. The perturbations also introduce 
ringing into the power spectrum which is retained in the 
full model. Because neighboring redshift modes produce 
similar responses in the power spectrum, there exists a 
degeneracy in the recovery of 5x(zi). 

Under the linear approximation, we can pose the in- 
verse problem as the recovery of 5x(z{) from the linear 
response 



scf E = c, 
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TiiSx(zi) + N e . 



(2) 



FIG. 2: Top: _E-mode polarization power spectrum for: the 
fiducial model of Fig.0(thick); the step function model (thin); 
the step function model with deviations transferred onto the 
fiducial model (dashed); instrumental noise Wp 1 ^ 2 (denoted in 
/iK-arcmin) that roughly brackets expectations from WMAP 
and Planck (long dashed). Bottom: the transfer function or 
fractional power spectrum response to a delta function per- 
turbation of unit amplitude at 8 < z% < 25. 



where Ng represents noise sources. Instrumental noise 
levels that roughly bracket the expectations for WMAP 
and Planck [l2| are shown in Fig. but systematic ef- 
fects such as residual foreground contamination will likely 
dominate the errors. With a characterization of the sta- 
tistical properties of the noise and a regularization of the 
reconstructed signal, e.g. by placing a smoothness cri- 
teria on x(zi) through a two-point prior, standard tech- 
niques such as Wiener filtering allow for model recon- 
struction. It is important to take a fiducial model that 
is close to the observed power spectrum since linear re- 
sponse will only hold for small perturbations. It will also 
be important below in that the cosmic variance of the 
assumed signal sets the fundamental noise in the recon- 
struction. We chose a fiducial model with an ionization 
x(zi) = 0.5 in the parameterized regime so that, neglect- 
ing helium reionization |6J, the maximal excursion from 
the model is Sx (zi) = 0.5. 

As a worst case scenario consider the reconstruction 
of the step function ionization in Fig. ^ from the fidu- 
cial model, Sx(zi) and the transfer matrix. Here all of 
the parameters x(zi) are offset from the fiducial choice 
by the maximum value. The predicted power spectrum 
is shown in dashed lines in Fig. [21 (top). Note the agree- 
ment is still quite good at the peak of the power but 
fails in the low signal regime where small errors in the 
transfer function lead to large fractional effects. Still, 
given realistic measurement noise and an iterated choice 
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FIG. 3: Eigenmodes. (a) 5 best (decreasing eigenvalue thick 
to thin) constrained eigenmodes or linear combinations of ion- 
ization history, (b) 5 worst constrained eigenmodes. 

of fiducial model, the linear response approximation can 
provide useful tools for representing the data. A detailed 
consideration is beyond the scope of this work. Instead 
we will show that a modification of the forward approach 
suffices to extract essentially all of the information in the 
CMB. 



III. CMB OBSERVABLES 

To better quantify the information contained in the 
power spectrum, let us consider the ultimate limit of an 
all-sky experiment that is cosmic variance limited. The 
variance of the power spectrum is then given by 

{8Cf E 8Cr) = ^(Cn%« (3) 

and hence the covariance of the ionization parameters 
(Sx(zi)Sx(zj)) ps (F _1 )jj, where 

Fij =$> + l/2)T«T^ (4) 

e 

is the Fisher matrix. The structure of the transfer matrix 
implies a large covariance between estimates of 5x(zi) 
and renders the delta-function representation difficult to 
visualize. 

Consider instead the principal component representa- 
tion based on the orthonormal eigenvectors of the Fisher 
matrix, decomposed as 

For a fixed the Sip specify linear combinations of the 
Sx(zi) for a new representation of the data 

m p = 2J S ifi Sxi , (6) 
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FIG. 4: Eigenmode statistics. Top curve: rms error cr M on 
mode amplitude; dashed line represents a physicality prior on 
x; only the first 5 modes contain interesting information. Bot- 
tom curve: optical depth per unit-amplitude mode r M . Middle 
curve: rms error on total optical depth shown as the cumula- 
tive contribution from modes < fi; dashed line represents the 
physicality prior on x. 

where the covariance matrix of the mode amplitudes is 
given by 

(m^rriv) = a^S^ . (7) 

In other words, the eigenvectors form a new basis that 
is complete and yields uncorrelated measurements with 
variance given by the inverse eigenvalue. The largest 
eigenvalues correspond to the minimum variance direc- 
tions and the first 5 are shown in the upper panel of 
Fig. |3 The first two correspond essentially to the aver- 
age ionization at high redshift and low redshift respec- 
tively. The lower panel shows the directions with the 5 
highest variances. Here neighboring delta modes with 
similar responses compensate each other to leave the ob- 
servable power spectrum unchanged. The rms of each 
mode is shown in Fig. Because the ionization fraction 
cannot be negative and the amplitude of each mode is 
~ 0.5, only the first 5 modes with < 1 have useful in- 
formation. An added benefit of the principal component 
representation is that the structure in the lowest modes 
is invariant under refinement of the binning scheme Az. 
In £ space, the first mode controls the high i power, the 
second the low I power and the third through fifth adjust 
the ringing in the spectrum. 

These eigenmodes provide a good meeting ground be- 
tween observations and models. The amplitude of 
these 5 best modes may be added to the usual CMB pa- 
rameter estimation chain and the results compared to 
model predictions for m M without significant loss of in- 
formation. As an example, in Figure [S] we have repre- 
sented a complex ionization history (inset, thick-dashed) 
through its first 1 through 5 eigenmodes. For this ioniza- 
tion history the first 3 eigenmodes suffice to recover the 
observable power spectrum. The temperature polariza- 
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FIG. 5: Representation of an arbitrary ionization history with 
the first 5 eigenmodes. Inset: ionization history in the true 
model (thick dashed) compared with representation with 1 
to 5 eigenmodes (solid, increasing thickness) away from the 
fiducial model (thin dashed). Main panel: resulting predic- 
tion for the power spectrum. With three or more modes the 
prediction is indistinguishable from the true model. 

tion cross spectra converges similarly. Note that in the 
eigenmode analysis the assumed ionization fraction can 
go unphysically negative but this causes no difficulty for 
a Boltzmann code and no ambiguity once interpreted as 
a mode amplitude. 

Finally it is interesting to consider the implications 
for the total optical depth or equivalently errors on the 
initial amplitude of fluctuations a\ n g, rs <j t . A precise 
normalization of the initial conditions is crucial to dark 
energy studies utilizing the evolution of structure [l3| . 
Each principal component mode perturbs the total opti- 
cal depth by 

^=E^^y- ( 8 ) 

This quantity is shown in Fig. 0] (lower curve). The best 



constrained mode bears a r uncertainty of 0.0026. Note 
that although the rms uncertainty in the higher modes 
is large, the optical depth contributed by them is com- 
pensatingly small. Given that the total optical depth 
variation is St = t m to m , the total variance is given by 

u\ = ^2 T ^, a fi- We plot in Fig. 2] the cumulative variance 
from increasing the number of modes included in the sum. 
The contributions are again dominated by the first few 
modes and placing a prior of ay, < 1 against unphysical 
values of x is enough to eliminate the small contribution 
from the higher modes. Thus, even with an arbitrary 
ionization history, the total optical depth can in princi- 
ple be measured to o~ T ~ 0.01 by an ideal experiment. 
Note that this is an uncertainty due to cosmic variance 
and not a bias. In the example of Fig.|3 the total optical 
depth with 0, 2 and 4 modes is r = 0.1732, 0.1498, 0.1379 
compared with a true optical depth of r = 0.1375. 

IV. DISCUSSION 

We have studied the effect of an arbitrary ionization 
history on the CMB polarization power spectrum and 
shown that the information lies in 5 broadly distributed 
modes in redshift. We have taken a complete basis for 
a range of redshifts 7 < z < 26 but this can be ex- 
tended as needed. The amplitude of these modes repre- 
sent a compact description of the observable properties of 
rcionization in the CMB and can be easily predicted as a 
filtered version of any given model for the ionization frac- 
tion. They can therefore serve as a model-independent 
tool for data analysis and model testing. They also can 
ultimately remove any bias in the measurement of the to- 
tal optical depth, or equivalently the initial amplitude of 
fluctuations, leaving a residual uncertainty from cosmic 
variance at the 1% level. 
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